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Abstract 

The birefringence effect in the nuclear pseudoelectric field of matter and an external electric 
field for a particle (deuteron, nucleus) moving in a storage ring is discussed. The influence of the 
birefringence effect on the EDM measurement experiments is considered. The attention is drawn 
to the possibility to measure the spin-dependent amplitude of the elastic coherent scattering of a 
deuteron by a nucleus, the electric polarizability of a deuteron (nucleus). Using a gas target with 
polarized nuclei also allows to study P-,T-odd interactions. 

1 INTRODUCTION 

The phenomena of spin rotation and spin dichroism (birefringence effect) for particles with the spin 
S > 1 in an unpolarized medium were theoretically described for the first time in P3E]- The peculiarity 
of this phenomenon is in conversion of the vector polarization to the tensor one and vice versa (similar 
the linear to circular polarization conversion for a photon in an optically anisotropic medium - the well 
known optical birefringence). However, unlike the optical birefringence, the birefringence effect for 
particles appear in isotropic matter (and even the spin of matter nuclei is either zero or unpolarized 
!). Anisotropy is provided by the particle itself (a particle with the spin S > 1 and mass M^0 has 
the intrinsic anisotropy). Deuteron spin dichroism was recently observed with the 20 MeV accelerator 
0. Further investigations of this phenomenon are planned to be carry out with a storage ring and 
an external beam [3J . Observation of particle spin rotation and spin dichroism (the birefringence 
effect) with a storage ring requires reducing of (g — 2) precession frequency (g is the gyromagnetic 
ratio). This precession appears due to interaction of the particle magnetic moment with an external 
electromagnetic field. The requirement for {g — 2) precession cancellation also arises when searching 
for a deuteron electric dipole moment (EDM) by the deuteron spin precession in an electric field in 
a storage ringjSlEj. Some extremely interesting ideas providing noticeable cancellation for (g — 2) 
precession are proposed in jSJ [7] 

According to 6 balancing the energy of the particle and the strength of the electric field in the 
storage ring provides to reduce and even zeroize the (g — 2) precession frequency. As a result, EDM- 
caused spin rotation grows linearly with time jHIIZj- Note that, when (g — 2) precession is suppressed, 
the angle of spin rotation induced by the birefringence effect grows linearly with time, too. 

The effect of deuteron (nucleus) birefringence in matter reveals itself in a storage ring due to 
presence of the residual gas inside the storage ring and use of a gas jet (gas target) for deuteron 
(nucleus) polarization analysis. Moreover, birefringence also occurs in the solid target, which is used 
for analysis of polarization of the deuteron (nucleus) beam outside the storage ring. Therefore, the 



phenomenon of birefringence in the gas medium and polarimeter would appear as a systematic error in 
the EDM measurements j7]. In addition, study of the birefringence phenomenon is of self-importance 
since it makes possible to measure the spin-depended part of the forward scattering amplitude. 

Lastly, the action of the electric field on the deuteron rouses one more mechanism of deuteron 
spin rotation and oscillations (the phenomenon of birefringence in an electric field) conditioned by the 
deuteron tensor electric polarizability jS). 

In this paper several our preceding works concerning the birefringence effect in the nuclear pseu- 
doelectric field of matter and an external electric field ^ El El lU El El H51 are summarized for 
a particle (deuteron, nucleus) moving in a storage ring . The attention is drawn to the additional 
possibility to measure the spin-dependent amplitude of the elastic coherent scattering of a deuteron 
by a nucleus, the electric polarizability of a deuteron (nucleus). Use of a gas target with polarized 
nuclei also allows to study P-,T-odd interactions. 



2 The index of refraction and effective potential energy of particles 
in medium. 



Close connection between the coherent elastic scattering amplitude at zero angle /(0) and the refraction 
index of medium N has been established as a result of numerous studies (see, for example, [HUTU]): 



(1) 



where p is the number of particles per cm 3 , k is the wave number of a particle incident on a target. 

The expression (^Q) was derived in assumption that N — 1 < 1. If k — ► then (N — 1) grows and 
expression for N has the form 

iV 2 = l + ^/(0) 
Let us consider particle refraction on the vacuum-medium boundary. 

The wave number of the particle in vacuum is denoted k, k' = kN is the wave number of the 
particle in the medium. 
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Figure 1: Kinetic energy of a particle in vacuum is not equal to that in medium. 

As it can be seen, the kinetic energy of a particle in vacuum E = 4^- is not equal to that in the 
medium E' = h2 ^ N2 . 

2m 

From the energy conservation condition we immediately obtain the necessity to suppose that the 
particle in the medium possesses an effective potential energy (see the detailed theory in JU].) This 
energy can be easily found from the evident equality 

E = E' + V 



i.e. 



V = E - Ef 



27Th 2 



p/(0). 



(2) 



Above we considered the rest target. But in storage rings moving bunches can be used as a target. 
Therefore we should generalize the expressions ()1I2|) for this case. Thus, let us consider the collision 
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of two bunches of particles. Suppose that in the rest frame of the storage ring the particles of the first 
beam have the energy E\ and Lorentz- factor 71, whereas particles of the second beam are characterized 
by the energy E2 and Lorentz-factor 72. Let us recollect that the phase of a wave in a medium is 
Lorentz-invariant. Therefore, we can find it by the following way. Let us choose the reference frame, 
where the second beam rests. As in this frame particles of the second beam rest, then the refraction 
index can be expressed in the conventional form l[T|): 

N[ = l + ^f(E>,0), (3) 
k[ 

where p' 2 = 7 2 _1 /°2 is the density of the bunch 2 in its rest frame and P2 is the density of the second 
bunch in the storage ring frame, k[ is the wavenumber of particles of the first bunch in the rest frame 
of the bunch 2. Let us denote the length of the bunch 2 in its rest frame as L, L = 72 /, where I is the 
length of this bunch in the storage ring frame. 

Now we can find the change of the phase of the wave caused by the interaction of the particle 1 
with the particles of bunch 2: 

<f> = k[(N[ - l)L = ^f(E[,0) L = ^/(^ j0 )A4 I, (4) 

It is known that the ratio ^ ^} is invariant, so we can write ^^}'°^ = , where f(E,0) is 

the amplitude of elastic coherent forward scattering of particle 1 by the moving particle 2 in the rest 
frame of the storage ring. 

As a result 

= ^/(^,0) I = *¥*f(E lt 0) v rel t, (5) 
k\ k\ 

where v re i is the velocity of relative motion of the particle 1 and bunch 2, if the first particle is 
nonrelativistic, whereas the second one is relativistic, then v re i ~ c, where c is the speed of light, t is 
the time of interaction of the particle 1 with the bunch 2 in the rest frame of the storage ring. 

The particle with the velocity v\ = ^ passes the distance z = v\t over the time t. It should be 
noted that the path length z differs from the length of bunch 2, because it moves. Expression © can 
be rewritten as: 

0=^/(^,0) '^z = k 1 (N 1 -l)z, (6) 
k\ v\ 

where the index of refraction of the particle 1 by the beam of moving particles 2 is: 

N, = l + ^^/(E,0) (7) 
ki vi 

If V2 = 0, the the conventional expression (1) follows from (JT)) 

So, let us consider particles captured to a trap. They are nonrelativistic. From (J2J) the effective 
potential energy V can be expressed as: 

27rfi 2 r 

V = -—p 2 -f(E 1 ,0)= (8) 
= —2-Khp2C - 
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where 72 is the Lorentz-factor of the bunch 2 and the Lorentz-factor of the particle 1 is 71 = 1 
hk[ = TO1C72. 

It should be noted that the amplitude of coherent elastic scattering at zero angle depends on the 
T-matrix as follows: 

f{E [, 0) = -MMr(^) = (2n) 2 m l72 T(E[) (9) 
V1 + V2 

i.e. the amplitude of forward scattering is proportional to the Lorentz-factor 72. As a result, the 
quantity j^ 1 f(E[,0) depends on the energy of the particle only due to possible dependence of T- 
matrix on energy. From ([9110)1 one can obtain 

V = (2^fp 2 T{E[) (10) 



3 THE PHENOMENON OF ROTATION AND OSCILLATION OF 
DEUTERON (NUCLEAR) SPIN IN UNPOLARIZED MATTER 
(BIREFRINGENCE AND SPIN DICHROISM) 

Let us consider the refraction of a particle with the spin S > 1 in matter. 

According to 012] the index of refraction for such a particle depends on the particle spin and can 
be written as: 

^ = 1 + fr/(°)' ( n ) 

where/ (0) = TrpjF(0), p is the density of matter (the number of scatterers in 1 cm 3 ), k is the 
deuteron wave number, pj is the spin density matrix of the scatterers, F (0) is the operator of the 
forward scattering amplitude, acting in the combined spin space of the particle and scatterer spin J. 

For a particle with the spin 5=1 (for example, deuteron) in an unpolarized target /(0) can be 
written as: 

/(0) = d + di{Snf. (12) 

where S is the deuteron spin operator, n is the unit vector along the deuteron momentum k. 

The axis of quantization z is directed along the particle wave vector k. Considering only strong 
interactions, which are invariant to the parity transformation and time reversal, we may omit the 
terms containing S in odd degrees. Therefore, the refractive index for deuterons 

N =l + ^(d + d 1 S 2 z ) (13) 

depends on the deuteron spin orientation relative to the deuteron momentum. 

The refractive index for a particle in the state, which is the eigenstate of the operator S z of spin 
projection on the axis z: 

N = 1 + (d + d im 2 ) , (14) 

m is the magnetic quantum number. 

According to Eq.(|14j). the refractive indices for the states with m = +1 and m = —1 are the 
same, while those for m = ±1 and m = are different (3ftiV(±l) / 5ftiV(0) and Q?iV(±l) ^ 3?7V(0)). 
This can be obviously explained as follows (see FigEJ): the shape of a deuteron in the ground state is 
non-spherical. Therefore, the scattering cross-section ctm for a deuteron with m = ±1 (deuteron spin 

is parallel to its momentum k) differs from the scattering cross-section a± for a deuteron with m = 0: 

°\\ ^ cr ± 9tfj,(0) = ^<7|| + 9f/ ± (0) = ^a ± . (15) 
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Figure 2: Two possible orientation of vectors S and n = £: a) S\\n; b) 5_Ln 



According to the dispersion relation K/(0) ~ $(9/(0)), then 5ft/j_(0) ^ 5R/j|(0). 

From the above it follows that deuteron spin dichroism appears even when the deuteron passes 
through an unpolarized target: due to different absorption, the initially unpolarized beam acquires 
polarization or, yet more precisely, alignment. 

Let us consider the behavior of the deuteron spin in a target. The spin state of the deuteron 
is described by its vector and tensor polarization (p = (S) and p^ = (Qik), respectively, here 

Qij = | (^SiSj + SjSi — §<%Y When the deuteron moves in matter its vector and tensor polar- 
ization appears changed. To calculate p and pik one need to know the explicit form of the deuteron 
spin wave function ip 

The wave function of the deuteron that has passed the distance z inside the target is: 

V> (z) = exp (^ikNz^j ipo, (16) 

where ipo is the wave function of the deuteron before entering the target. The wave function tp can 
be expressed as a superposition of the basic spin functions Xm> which are the eigenfunctions of the 
operators S 2 and S z (S z Xm = mXm)'- 



Therefore, 



m=±l,0 



a l \ / ae iS le ikN 1Z 

q, = [ a | = | be i5 °e ikN ° z 

ce iS-i e ikN~iz 



ae iSi e ikNiz 



(18) 



^ e i8 0e ikN z 
ce iS-i e ikNiz 



(according to the above Ni = iV_i). 

Suppose plane (yz) coincides with the plane formed by the initial deuteron vector polarization 
Po ^ and the momentum k of the deuteron. In this case 5i — 5q = 5q — <5_i = ^, and the components 
of the polarization vector at z = are p x = 0,p y ^ 0, and p z ^ 0. 

-> (ty\s\ty) 

The components of the vector polarization p = (S) = inside the target are: 



5 



Py 



v/2e-^( ff °+ CT iH( a - c) sin 
y2 e -|p(<ro+ai) Z6 ( a + c ) cog fesft^ 



(* I *) 



(a 2 - c 2 ) 



(* I *) 



(19) 



Similarly, the components of the tensor polarization Qij = | ( SiSj + SjSi — ^5ij ) are expressed as: 



Pxx 

Pyy 

Pzz 
Pxy 

Pxz 
Pyz 



■\ (a 2 + c 2 ) e- paiz + b 2 e- paoZ - 3ace- pai " 
(* I *> 

\ (a 2 + c 2 ) e- paiZ + b 2 e- pCT0Z + 3ace- pai " 

(* | *) 
2 2 + c 2 ) e~ paiz - 2b 2 e-P aoz 



0, 



_^ e -Ap(«o+<n)* b ( a + c ) sin C^fiKdi 
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V2 



(20) 



where (* | *) = (a 2 + c 2 ) e"^ 2 + b 2 e~ p ^ z , a = f 3f/ , ^ = £ / = d, ft = d + d x . 

According to ()19[). (|2U|) spin rotation occurs when the angle between the polarization vector p and 
momentum k of the particle differs from ^. 

For example, when 3tdi > 0, the angle between the polarization vector and momentum is acute 
(see Fig|3J) and the spin rotates left-hand around the momentum direction, whereas an obtuse angle 
between the polarization vector and the momentum gives rise to right-hand spin rotation (see FigEJ. 
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Figure 3: 
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When the polarization vector and momentum are perpendicular (a transversely polarized particle), 
then the components of the vector polarization at z = are: p x = 0, p y ^ 0, and p z = 0. In this case 
a = c and dependance of the vector polarization on z can be expressed as: 

Pi 



Py 

Pz 
Pxx 

Pyy 



= o, 






^p(°o+°i)z 2bacos hn> miZ \ 




<* 1 *> 


= o, 




-4a 2 


e —(xrt» _|_ ^Sg-pcro^ 




(* 1 *> 


2aV 


_|_ y2 e -pa- z 




(* | *) 


2aV 





(21) 



^ e -|p( CT o+^)^2a6sin (^^z 



0. 



According to ()22|) the vector and tensor polarization oscillate when a transversely polarized deuteron 
passes through matter (see FiglSJ) i.e. the vector polarization converts to the tensor one and vice versa 
(similar the linear to circular polarization conversion for a photon in an optically anisotropic medium 
- the well known optical birefringence). However, unlike the optical birefringence, the birefringence 
effect for particles appear in isotropic matter (and even the spin of matter nuclei is either zero or 
unpolarized !). Anisotropy is provided by the particle itself (a particle with the spin S > 1 and mass 
M ^ has the intrinsic anisotropy). 
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Figure 5: 



According to the above (see (j2J) a particle possesses some effective potential energy V in matter. 
If the particle spin S > 1, then this energy depends on the spin orientation 12 Q] 

^— ^/f"). < 22 > 

where M is the particle mass, /(0) is the spin dependent zero-angle elastic coherent scattering am- 
plitude of the particle, iV is the density of the scatterers in the medium (the number of scatterers in 
lcm 3 ), 7 is the Lorentz factor. Substituting /(0) for a particle with the spin S = 1 in 1)22(1 in the 
explicit form one can obtain El 13 

V = -^-N (d + d^SnV) , (23) 



where n is the unit vector along the particle momentum direction 
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Let the quantization axis z is directed along n and m denotes the magnetic quantum number. 
Then, for a particle in a state that is an eigenstate of the operator S z of spin projection onto the 
z-axis, the effective potential energy can be written as: 

V = -^N(d + d im 2 ). (24) 

According to (|24[) splitting of the deuteron energy levels in matter is similar to splitting of atom 
energy levels in an electric field aroused from the quadratic Stark effect. Therefore, the above effect 
could be considered as caused by splitting of the spin levels of the particle in the pseudoelectric nuclear 
field of matter. 

Let an electric field acts on a deuteron (nucleus). The energy Ve of the deuteron beam in the 
external electric field E due to the tensor electric polarizability can be written in the form 

V E = -^a ik EiE k , (25) 

where ctik is the deuteron tensor electric polarizability, E{ are the components of the electric field. 
This expression can be rewritten as follows: 

V E = a s E 2 - a T E 2 {Sn E ^ ^ , (26) 

where as is the deuteron scalar electric polarizability, ay is the deuteron tensor electric polarizability, 
He is the unit vector along E. 

Comparison of (|26|) with 1)23(1 provides to conclude that in an electric field we can observe the 
effect of spin rotation and oscillations about the E direction for a particle with S > 1 , too [S]. 



4 THE EQUATIONS FOR THE POLARIZATION VECTOR AND 
QUADRUPOLARIZATION TENSOR OF THE DEUTERON 
BEAM IN A STORAGE RING 

Considering evolution of the spin of a particle in a storage ring one should take into account several 
interactions: 

1. interactions of the magnetic and electric dipole moments with an electromagnetic field; 

2. interaction l|2fi|) of the particle with the electric field due to the tensor electric polarizability 

3. interaction (|23|) of the particle with the pseudoelectric nuclear field of matter. 
Therefore, the equation for the particle spin wavefunction is: 

where *&(t) is the particle spin wavefunction, Hq is the Hamiltonian describing the spin behavior 
caused by interaction of the magnetic moment with the electromagnetic field (equation (|27|) with the 
only Hq summand converts to the Bargman-Myshel-Telegdy equation), describes interaction of the 
deuteron (nuclear) EDM with the electric field. 

Let us describe motion of a deuteron in a storage ring in external magnetic and electric fields. 
Particle spin precession induced by interaction of the magnetic moment of a particle with an external 
electromagnetic field can be described by the Bargman-Myshel-Telegdy equation 

f=[pxfU (28) 
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where t is time in the laboratory system, 



O = — 

mc 



a + — ] B 

7, 



7+1 



P-B)0 



7 + 1 



x E 



(29) 



m is the mass of the particle, e is its charge, pis the spin polarization vector, 7 is the Lorentz- factor, 
P = v/c, v is the particle velocity, a = (g — 2)/2, g is the gyromagnetic ratio, and if are the electric 
and magnetic fields in the point of particle location. 

If a particle possesses an intrinsic dipole moment then the additional term that describes the spin 
rotation induced by the EDM should be added to ((28|) 6 



dp, 



edrn 



dt 



d 

h L J 



px 1/3 x B + E 



(30) 



where d is the electric dipole moment of a particle. 

As a result, evolution of the deuteron spin due to the magnetic and electric momenta can be 
described by the following equation: 



dp 
~dt 



e 








px | 


R) 


mc 





7 



7 + 1 



P'B)P 



7 



7 + 1 



(3x E 



+ d 



p x 



x B + E 



.(31) 



According to the section 2, the equation (|3*Tj) does not describe the particle spin evolution in a 
storage ring completely. The expression (|3*Tj) should be supplemented with the additions given by 
interactions Ve and V (see ([23127)) ) . 

This additional contribution could be found by the aids of the particle spin wavefunction ^>{t) (see 
EZJ)). 

The equations describing the time evolution of the spin and tensor polarization caused by the 
phenomena of birefringence can be written as: 



dp 

~dt 

dpik 
dt 



d {^(t)\S\*(t)) 

dt {mm)\ ' 

d {9(t)\Q ik \*{t)) 

dt m)Mt)\ 



(32) 



(S{Sk + SkSi — 3<5jfci) is the tensor of rank two 



where ^(t) is the deuteron wave function, ■ 
(the tensor polarization). 

The equations (|32|> contain initial phases that determine the deuteron wave function. Therefore, a 
partly polarized beam can not be described by such equations. So the spin density matrix formalism 
should be used to derive equations describing the evolution of the deuteron spin [TB] . 

The density matrix of the system " deuteron+target" is 



P = Pd® Pt, 

where pd is the density matrix of the deuteron beam 



Pd = I(k) 



-I + ^p(k)S + ^p ik (k)Qik 



(33) 



(34) 



I(k) is the intensity of the beam, p is the polarization vector, pik is the tensor polarization of the 
deuteron beam, pt is the density matrix of the target. For an unpolarized target pt = I, where I is 
the unit matrix in the spin space of target particle. 

The equation for the deuteron beam density matrix can be written as: 



dpd 
dt 



H,p d 



+ 



dpd 
dt 



(35) 



col 
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where H = H + V d + V E , 



Vd 
Ve 

He 



d(j3xB + E)S, 



(36) 



_A 2 



The collision term ( 



col 



= a s i/3 x B + Ej — (xt y(3 X B + E) lSn E ) , 

E + (3 x B 
\E + 0xB\' 

can be found by the method described in |12j : 
^ [F(0 = 0)p - P F + {0 = 0)] + J dnF(k')p(k')F + (k') 



(37) 



where = k + q, qis the momentum carried over a nucleus of the matter from the incident particle, 
v is the speed of the incident particles, N is the atom density in the matter, F is the scattering 
amplitude depending on the spin operators of the deuteron and the matter nucleus (atom), F + is 
the Hermitian conjugate of the operator F. The first term in (|37jl describes coherent scattering of a 
particle by matter nuclei, while the second term is for multiple scattering. 
Let us consider the first term in (|37|): 



(i) 



col 



vN- 



,2m 
~~k~ 



f(o) Pd - Pd f(py 



(38) 



The amplitude /(0) of deuteron scattering in an unpolarized target at the zero angle is 

/(0) = S Pt F(0)p t . (39) 

This amplitude can de rewritten according to ()23l) as 

/(0) = d + d 1 (Sn) 2 , (40) 

where n = k/k, k is the deuteron momentum. 
As a result one can obtain: 



Finally, the expression (|35|) reads 



dpd 
dt 



^[VPd- PdV + ) + vNS Pt I dnF(k )p(k )F + {k 



(41) 



(42) 



The last term in the above formula, which is proportional to Sp t , describes the multiple scattering 
process and spin depolarization aroused from it. Henceforward we consider such time of experiment 
(such effective length for a particle in a matter) that provides to neglect this term. 
The intensity of the beam is 



SPdPd- 



Consequently 



dl 



2-Ki 



vN—S Pd [f(0)p d - Pd f + (0)] 



(43) 



(44) 
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Substituting J3U) and (jHH) into (jUJ) we can get 

dl X 



dt 



2+p ik nin k ] I(t) + al(t), 



(45) 



where x 



AttvN 



Imdi 



vN(a± — o"o), a 



AirvN 



-vNao, o\ and ctq are the total cross- 



sections of deuteron scattering by a nonpolarized nucleus for the magnetic quantum numbers m = 1 
and m = 0, respectively. 

The vector polarization of the deuteron beam p is determined as 



P 



SpdPdS Sp d p d S 



SpdPd I 

From l)46|) one can get the differential equation for the beam polarization 

dp _ Sp d (dp d /dt)S _ _ ^Sp d (dp d /dt) 
dt ~ I{t) P I(t) 

The expression for the tensor polarization is 

Sp d PdQik Sp d p d Qik 



Pik 



SPdPd 



I 



(46) 



(47) 



(48) 



where Q ik = § (SiS k + S k Si - ^5 ik i). 

The change of the tensor polarization can be written as 



dp 



ik 



Sp d (dp d /dt)Q ik Sp d (dp d /dt) 
Pik- 



(49) 



dt I(t) I(t) 

Using (|34[) and (|2<S|) . (j47j) and (|49|) we can get the equation system for the time evolution 



of the deuteron polarization vector and quadrupolarization tensor (n = k/k, n E 

Pxx + Pyy + Pzz = 0) 16 



E+f3xB 
\E+P*B\' 



p X 



{ (« + B ~ (0- B) - (f - ^) /3 x 



dp 
dt ~ 

+i 



+ |[n x n] — ^fp— \ (u ■ n)p 



+ 



pxlE + f3xB) +%(n(n-p) +p)+ 



z f 3 v •« ji> 3 n [nE><n E \ 



^ d l — {^jkrPij^r + EjirPkj^r) + 



(50) 



{~5 + nifl k + \Pik ~ \{™i n k + nin k ) + \{n- n)5 ik } + 
+¥ ([" x pW + n;[n x p] fc ) - l(n • n)p ik - 



+X ([" x + "<i[n x p] fc ) - |(n • n )p. 
'\ aT h U ^ E x P\i n E,k + n E ,i[n E x p] fc ) , 

where -E e // = (E + f3 x B), rj = —^j^-Redi, n i = pi k n k , n E i = pi k n Ejk , f2 r (d) are the components of 
the vector Q(d) (r = 1, 2, 3 corresponds x, y, z): 



Q{d) = — 
mc 



a + - I B - a- 

7 



7 



7 + 1 



P'B)P 



7 + 1 



x E\ + 



(51) 



Then we consider the spin rotation about the particle momentum. According to Ej the spin 
precession caused by the magnetic moment {{g — 2) precession) can be minimized and even zeroized 
by applying a radial electric field. 
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The angles of spin rotation caused by both the EDM and birefringence effect are small for the 
considered experiment duration. Therefore, the perturbation theory can be used for (jST)|) solution. 



P(t) =p° + — 
mc 



p° x {aB + 



1 



a \(3 x E 



d 



p u x ( E + (3 x B 
2x 



7 2 -l 

t + {|(n(n • p°) +p°)t + |[n x n' ]t - 



t + 



-fi?t - |(n • n )P°*} " l^^f^ x n' E0 ]t, 



(52) 



= Pik ~ ^{tjkrPij + ZjirPkj) \ aB + 



a ) (3 x E \ t 



(£jkrPij + ZjirPkj) [E + x B\ t + 



+{x 



1 



1 



,0 



1 



1 



t + 



+ n ink + -p" ik - -{n i0 n k + rnn k0 ) + -(n • n )8 ik 

+ -j {[n x p\n k + rn[nx p%) t-~(n- n )p? k t} - 
3 axEl ff 



(53) 



where p° is the beam polarization at to = 0, n i0 = p® k n k , n E0 i = p® k nE,k, Pi k are the components of 
the tensor polarization at the initial moment of time. 

In real situation, even when {g — 2) precession is suppressed, nevertheless, the rotation angle can 
appear large enough (during the experiment the spin can rotate several turns [7]). Absorption can 
also appear significant. In this case one should analyze the system (|50|) instead of perturbation theory 
results (J52|l. 

Thus, according to (f52*|) . (|53*)) the spin behavior of a deuteron rotating in a storage ring is caused 
by several contributions: 

1. spin rotation which is described by Bargman-Myshel-Telegdy equation; 

2. rotation due to the deuteron EDM; 

3. rotation and dichroism due to the the phenomena of birefringence in matter; 
and 

4. spin rotation due to the phenomena of birefringence in an electric field. 

Let us consider some particular cases. 

Case I. Suppose the vector polarization is parallel to the the z axis, i.e. p® = p® = 0, p® ^ 0, 



Pik 



0, if % + k, p xx + 0, p£, + 0, pi 



0. 



p 



Figure 6: The original orientation of the polarization vector 
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Suppose {g — 2) spin precession caused by the magnetic moment is zero, then one can obtain 



Py{t) 



Px(t) = 
p°(E + f3B)t, 



Pz(t) 



Pxx(t) ■ 
Pyy{t) 



fxx 
,0 



+ 



X 
3 

X 



Pz + \xPzt, 



Pyy 3 



Pzz(t) 



l +Pyy) t ' 
2 X, 



(54) 



d 



Pxz(t)=0, 



Pyz(t) = TPyy (E + f3B)t, 



n 

Pxy 



3 OL T Ef 



e //„o 

I z 



The solution (|55|) shows that even in this practically ideal case (when the polarization vector is exactly 
parallel to n) change of polarization due to the birefringence effect leads to the appearance of additional 
components of p and pn~ along with the components aroused from the deuteron EDM. 
Case II. Let us consider now the more real case. 

Suppose the angle between the initial polarization vector and the z axis is acute (FigEJ) Let us 



Figure 7: The original orientation of the polarization vector 



express the effect caused by the gas inside the storage ring. Then the solution of H52I53|) can be written 
as follows: 

u\ o X o , X o o 4. ^ o * 

Pxit) = p x - - p x t - - p zz Pxt-T^ Pyz* 

Py® =Py-l(E + pB)p«t-^p«t-^ p o p o t+ V_ p zf 

Pz(t)= P ° + ^(S + /3B)pft + |p»*-| p ° p» t 
Pxx (t) = pjL - | (1 + O t - | p° zz pl x t 
Pyy® = Pyy ~ \ Py* (E + f3B)-^{l+ p° xx ) t - | p° t p° yy t (55) 



Pzz(t) = P° zz + 1 d - p yz (E + (3B) + | (2 - pi) t - | p^J t 

Pxy^) = Ply - l PXZ {E + PB) + l Ply t - l Pi-, Ply t 

p xz (t) = pl z + | Pxy (E + (3B)-\ Plz t ~ y Pyt ~ § Pzz Plz t 



?yz + I (Pyy ~ Pzz) (E + ?B) - | p° z t +^p°t-^ p Q zz p° z ; 
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According to (j5fij) the change in components of vector and tensor polarization caused by the EDM 
is mixed with the contributions from the birefringence effect to the same components. In much the 
same mixing of contributions from the EDM and the tensor electric polarizability ay appear, but here 
we does not dwell on them because of the expressions unhandiness. 

Thus, the changes in the deuteron vector and tensor polarization are the result of several mecha- 
nisms: 

- rotation of the spin in the horizontal plane (E,n) (FigEJ)- 

B 



n 



E 



Figure 8: The respective orientation of E, B, and n 

This rotation is risen due to interaction of the magnetic moment with external fields. The rotation 
frequency is expressed as: 



w„ = — \aB + f -5— - a ) x E } ; 
mc I \7 Z — 1 



(56) 



rotation of the spin in the vertical plane (B,n) caused by the electric dipole moment (Fig|UJ); 

J 




Figure 9: Rotation of the polarization vector due to birefringence in an electric field 
The rotation frequency is 



(57) 



- rotation caused by the phenomenon of birefringence in a medium, this is precession in the vertical 
plane (B,E) (FigHHJ); 



B 



p 



n 



Figure 10: Rotation of the polarization vector due to the phenomena of birefringence 
The rotation frequency is 



2vriV 
M7 
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hRedi. 



(58) 



- rotation due to the phenomenon of birefringence in an electric field in the vertical plane (B, n) 
(Fig lllJI . i.e. in the same plane as the rotation caused by the EDM. 




Figure 11: Rotation of the polarization vector due to birefringence in an electric field 



The rotation frequency is 



CUE 



eff 



(59) 



Besides the above rotations the transitions from the vector polarization into the tensor 
one and spin dichroism occur. 

Moreover, the spin dichroism leads to the appearance of the tensor polarization. 
Let us compare the frequency and the angle of polarization vector rotation caused by the EDM 
with those caused by the birefringence effect. 

1. The spin rotation frequency caused by the EDM is determined by the formula 1|30|) 



dE d QJJ 

Uedm = — r~ + -rpB. 

n n 



(60) 



We can get co e d m ~ 3 • 10 7 rad/s for the storage ring with E = 3.5 MV/m, B = 0.2 T and expected 
value of the deuteron EDM d ~ 10 e ■ cm and oJedm ~ 

3 • I0~ 9 rad/s for EDM d ~ 10" 29 e • cm. 



2. The spin rotation frequency caused by the phenomena of birefringence in a residual gas: 



2ttN _ , 
u = ~r~z — n Recti. 



M7 



(61) 



Using the last formula one can get uj ~ 2 • 10~ 7 rad/s for N = 10 9 cm -3 (suppose the pressure inside 
the storage ring ~ 10~ 7 Torr), Redi ~ 10 -13 . This effect depends on the density N (depends on the 
pressure inside the storage ring). 

3. The spin rotation frequency caused by the phenomenon of birefringence in the gas jet (gas 
target), which is used for the beam extraction to the polarimeter [J] (Fig|12|): 

In this case the effective rotation frequency (or the rotation angle for r = Is) is determined as 
follows: 

2irN t 

Uteff 



ft 



k 



■ I Red± v, 



(62) 



where uJteff is the effective frequency of spin rotation , <pt is the deuteron spin rotation angle for 
r = Is, Nt is the density of the target, I is the length of the target, v is the frequency of beam rotation 
in the storage ring, k = is the deuteron wave number. 



Coulomb nuclear 
scattering scattering 




gas jet 



segmented detector 



Figure 12: Polarimeter [7j 
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Really, the frequency uJt of the spin rotation in matter is 



<*>t = ; r 1 hRedi (63) 

The spin rotation angle 6 T = cJt T ti where n = ~ is the deuteron flying time in the target. The angle 
of rotation at 1 second is 9 = utTtU (during a second the deuteron passes through the gas target v 
times). As a result, 

U te ff = 0J t T t V, (64) 

so one can get 

2irN t . , I 2vrj 
Uteff = ^ri—nRedi-v=— i —Rediv, (65) 
JJ M7 v k 

where j = Nfl. Then using the experimental parameters j = 10 15 cm -2 , v ~ 10 5 - 10 6 we have 
uteff » 1(T 5 - W~ 4 rad/s. 

So the angle of polarization vector rotation for t = Is caused by the phenomenon of birefringence 
in the gas jet of polarimeter appears by two orders of magnitude greater than the angle of rotation 
due to the EDM. The additional contribution to spin rotation is also provided by the solid carbon 
target of polarimeter (see Figure lT2|) . 

4. The estimation for the spin rotation frequency caused by the birefringence in an electric field 
according to Q57JI for ay ~ 10 -37 cm 3 and E = 3.5MV/m is loe ~ W~ e rad/s. 

Let us estimate the value of spin dichroism. This characteristic is given by the parameter 
X = — v N(a\ — do) (see expressions (|52|) . (|53|0: 

- in the case of the scattering by the residual gas we have |x| ~ 0.5 • 10~ 6 for iV ~ 10 9 cm~ 3 ; 

- if a deuteron passes through the gas jet (g get) for r = Is is 

Xt = iOl - 00) v, (66) 

then for j = 10 15 cm~ 2 and v ~ 6 • 10 5 one can get xt ~ 3.4 • 10~ 5 . So we can conclude that there is 
the significant beam spin dichroism. 

It should be especially mentioned that the polarimeter (Fig |12|) includes a solid target. Therefore, 
the birefringence effect in the matter of the solid target arises, too. 

The possible influence of P-odd spin rotation and spin dichroism (i.e. different absorbtion cross- 
sections in matter for a deuteron with the spin parallel and antiparallel to the momentum <7ft 7^ <7fj) 
aren't to be forgotten. 



5 CONCLUSION 

The above analysis shows that the phenomena of deuteron birefringence in matter and an electric 
field should be studied in experiments for the EDM search with a storage ring. Moreover, they could 
imitate the spin rotation due to the EDM. Study of the birefringence effects in such experiments could 
provide to measure both the spin-dependent part d\ of the amplitude of the coherent elastic scattering 
of a deuteron by a nucleus at the zero angle and the tensor electric polarizability of a deuteron. 

It should be also mentioned that if the nuclei in the gas jet are polarized, then according to 0] 
the P-,T-odd spin rotation and dichroism appear in the storage ring. They are caused by the T- 
odd nucleon-nucleon interaction of a deuteron with a polarized nucleus and, in particular, interaction 
described by Vp t T ~ S \pn x ft], where pn is the polarization vector of gas target. 

P-even T-odd spin rotation and dichroism of deuterons (nuclei) caused by the interaction either 
Vt ~ {S\pN x n])(Sn) or ~ Sppj[([S x n]J)(Jn)] also could be observed [1] (here J > 1 is the spin 
of the polarized target nuclei, pj is the spin matrix density of the target nuclei). 
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6 Appendix 



Let us consider the amplitude of deuteron scattering by a proton depending on spin orientation with 
respect to the deuteron momentum j2j. 

For fast deuterons the scattering amplitude can be found in the eikonal approximation |131 114j . 
According to [Hj, the amplitude of zero-angle coherent scattering in this approximation can be written 
as follows: 

/(0) = ^-f (e tXD P*) - l) d 2 %(7>)|Vr (67) 

where k is the deuteron wavenumber, b is the impact parameter, i.e. the distance between the 
deuteron and the proton centers of gravity in the normal plane; tp ( r ) is the wavefunction of the 
deuteron in the ground state; \(p ( r )| is the probability to find proton and neutron (in the deuteron) 
at the distance '•r apart. The phase shift due to the deuteron scattering by a proton is 

XD = -— V D (b,z,-?±)dz' (68) 



hv 



CO 



~~r± is the ~r component perpendicular to the momentum of the incident deuteron, v is the deuteron 
velocity. The phase shift xd = Xi + X2, where \i an d Xi are the phase shifts caused by proton-proton 
and neutron-proton interactions, respectively. 

For the polarized deuteron under consideration the probability | ip ( r ) | is different for the different 
spin states of the deuteron. Thus for states with the magnetic quantum number m = ±1, the 
probability is \(p±i (r)\ , whereas for m = 0, it is \cpo ( r )| . Owing to the additivity of the phase 
shifts, the equation (|BT|) can be rewritten as 

/((.) - / V, ( ; ' --r) +/ 2 ( b+^A +2i tl (~b -%f)t 2 (-b + %f)\\ v (r)\ 2 dHd i r 



7T 



(69) 



where 



g*Xl(2) 



2i 

Attention should be given to the fact that the latter expression is valid if one neglects spin depen- 
dence of nuclear forces between the colliding proton (neutron) and proton. When this dependence is 
taken into account, the phase shift xd is an operator acting in the spin space of colliding particles and 
in the general case the expansion (|69j) is not valid. However, to estimate the magnitude of the effect 
of deuteron spin oscillations, spin dependence of nuclear forces may be neglected (spin dependent 
contribution was considered in ^0])- Let us also omit terms caused by the Coulomb interaction. From 
(HU) it follows 

/(0) = A(0) + / 2 (0) + ^ J t x (~b - ^)t2 (y + %f} \<p (?±, z)\ 2 d 2 bd 2 r ± dz (70) 

where 

k [\ r ~?\j2 t m D 



TT 



/i( 2 )(o) = - / t m ( t )d*z = f p(n) (0) 



and f p (n) (0) is the amplitude of the proton (neutron)-proton zero-angle elastic coherent scattering. 
The expression (|7T7|) can be rewritten as 



2ik 



/(0) = / 1 (0) + / 2 (0) + — / ti(0*2(?) <P[S-Tf,z) d^d 2 r)dz (71) 



2 
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Then from (|71jl 

oh r — > /— > \ 2 

u /(o)=«/i(o)+«/ 2 (o)-— »y ti(e)* 2 (^) ^(e -v,*) (72) 



2 A- 

9 /(0) = 9f /i(0) + 9 / 2 (0) + +— 3? / ti( £ )* 2 (?) 

7T 



d 2 £ d 2 ?? dz 



The spin oscillation period is determined by the difference in the amplitudes 5ft f(m = ±1) and 
sft /(m = 0). From (JSHJ) it follows that 

^ rf i = Z* 1 ^ ) <2 (^) P±i ( £ ~ ^' 2 ) V±i ( £ ~ z ) ~ ( £ - ~V> Z ) ( £ - d 2 £ d 2 ?7 dz 



di = / ti( £ )t 2 (~rf) P±i ( £ - ) z ) V±i ( £ - ~»f , z ) - Vo" ( £ - ~V, z) ¥>o ( £ - V, z 



d 2 £ d 2 ?? 



It is well known that the characteristic radius of the deuteron is large comparing with the range of 
nuclear forces. For this reason, when integrating, the functions t\ and t 2 act on (p as 5-function. Then 



AU 

SR d 1 = 9^(0) / 2 (0) / [>+ (0, z) p ±1 (0, z) - ^+ (0, z) ^ (0, z)] dz (73) 

^ Jo 

4k f°° 
$Sdi = — »/i(0) / 2 (0) / (0, z) p ±1 (0, z) - (0, z) (0, z)] dz 

n Jo 

The magnitude of the spin oscillation effect is determined by the difference 

[<pti (0, z) <P±i (0, z) - (0, z) ipo (0, z)] 

i.e. by the difference in distributions of nucleon density in the deuteron for different deuteron spin 
orientations. The structure of the wavefunction tp±i is well known: 



1 f u(r) 1 W(r) a . 



(74) 



where u{r) is the deuteron radial wavefunction corresponding to the S-wave; W(r) is the radial function 

corresponding to the D-wave; the operator §12 = 6(Slt r ) 2 — 25 ; ~n r = -y\ S = + ^2) and 

°" 1(2) a ^ e the Pauli spin matrices describing proton(neutron) spin. 
Use of dHJ) yields 



^d 1 = - JL Q {/i(0)/ 2 (0)} D = --^9 {/ p (0)/ n (0)} D 
3 dl = ^{/i(°)/2(0)}D = ^5ft{/ p (0)/ n (0)}D 

poo / 1 u(r)TV(r) _ 1 W 2 (r) 
\7/2 r 2 4 r 3 

where <r is the total scattering cross-section, one can obtain 

3 



where D = ^-^ u(r> ^2 (r> — | n Jf 1 ) dr, r 2 = £ 2 + z 2 . Applying the optical theorem 9 / 



9 di 



(K/ p (0K + 5ft/ n (0K)D 
3k 

■Ink V ' (47r) 2 "^V ~ v ^ /P ^ ^ ~ (2^F CTpCTn ' ° 



5ft di = 



2vr 2 
6 



(75) 



4^ ^ 



(76) 
(77) 



where (Tp( n ) is the total cross-section of the proton-proton (neutron-proton) nuclear scattering. 
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Figure 13: Rotation angle and dichroism in hydrogen target at high energies 




Figure 14: Spin rotation angle and dichroism in hydrogen target at deuteron energies lower than 5 
Gev. Solid curve corresponds to the calculation with the spinless N — N amplitudes 
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